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We present the design of a new type of compact toroidal, or "halo" , ion trap. Such traps may be 
useful for mass spectrometry, studying small Coulomb cluster rings, quantum information applica- 
tions, or other quantum simulations where a ring topology is of interest. We present results from 
a Monte Carlo optimization of the trap design parameters using finite-element analysis simulations 
that minimizes higher-order anharmonic terms in the trapping pseudopotential, while maintaining 
complete control over ion placement at the pseudopotential node in 3D using static bias fields. 
These simulations are based on a practical electrode design using readily-available parts, yet can be 
easily scaled to any size trap with similar electrode spacings. We also derive the conditions for a 
crystal phase transition for two ions in the compact halo trap, the first non-trivial phase transition 
for Coulomb crystals in this geometry. 
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I. INTRODUCTION 



Recent interest in novel ion trap geometries stems from 
both the need for smaller, more compact radiofrequency 
(RF) Paul ion traps [THS] and traps that have novel com- 
mon modes of motion [BHH] . We present a new compact 
"halo" ion trap geometry based on toroidal ion traps [9]- 
[TT] that might satisfy both of these needs. This new 
geometry has the advantage of being much smaller (on 
the order of a few hundred microns in diameter) than 
previous toroidal traps. Additionally, this geometry is 
very open both optically and electronically: ions trapped 
around the circular RF node interact via the Coulomb 
force across the circle with very little shielding, similar 
to ions in a Penning- type trap, giving rise to strong ion- 
ion interactions across the trap [12]. To date, Coulomb 
crystals have been studied in a number of different RF 
trap geometries including a linear trap [T3HT5] , a spher- 
ical trap [16 , ions confined to a two-dimensional plane 
[T7] , ions in a large toroid trap with no ion interaction 
across the circle [18-20 , and ions in a Penning trap [21]. 
Our compact halo trap geometry is complementary to 
many of these previously studied systems. This geome- 
try may also be suited to simulating small chemical rings 
[22] , circular clusters of electrons on the surface of liquid 
helium [23], and novel Ising models [24] . 

We begin by defining by the key parameters and char- 
acteristics of the compact halo ion trap geometry. Our 
setup is based on readily available electrode parts, mak- 
ing this configuration straightforward and inexpensive to 
implement. We model both the RF and static poten- 
tials in the trap using finite element analysis; we present 
the results from a Monte Carlo optimization of the ad- 
justable physical trap parameters yielding a trapping RF 
potential that is hyperbolic near the center of the trap. 
In order to provide complete control of the ions, we also 
optimize offset static potentials that provide control over 
the trap aspect ratio. In the final section we derive the 
conditions for a crystal phase transition for two ions in 
this new geometry. We show that, unlike a linear trap, 
there is a non-trivial phase transition for only two ions. 



II. HALO TRAP GEOMETRY 

A traditional RF quadrupole mass spectrometer or ion 
trap is created from four long cylindrical rods placed at 
the corners of a square as in Fig. pTa). Alternating RF 
potentials, d=Vo> are applied to the electrodes forming a 
2D quadrupole potential in the center of the square with 
a node that runs the length of the rods. A toroidal RF 
trap is formed by wrapping the ends of linear trap around 
on themselves, forming four concentric loops stacked two 
on top of the others as in Fig. [lib) [9j [11]. We pro- 
pose one further modification that shrinks the toroidal 
trap and simplifies trap fabrication. Since the inner rings 
are equipotential surfaces both on the top and the bot- 
tom, we shrink these rings down to solid, cylindrical elec- 
trodes, shown schematically in Fig. []Jc). The outer two 
electrodes are also modified: we use conductive tubes to 
create these potentials, again concentric with the central 
electrodes. We envision an ion trap with an RF node 
diameter of less than 1 mm. 

In order to control the relative trap aspect ratio in the 
s — z plane defined in Fig. [2j key to studying the crystal 
phase transitions [T3UT5] , we utilize both RF and static 
trapping fields. Figure [2] shows a cross-section of the 
halo geometry: rotating the s — z plane about the axis 
located at the left of the figure yields the 3D trap con- 
figuration illustrated in Fig. [3] The electrode geometry 
described above forms a quadrupole potential like that 
illustrated in Fig. [2|a) for potentials ±Vq applied to the 
center electrode and the outer tube. The effective RF 
potential, or pseudopotential, is illustrated in Fig. |2|b), 
the center of which is located along the z = center line 
of the geometry. To lowest order, the pseudopotential 
is symmetric in the s — z plane; we discuss higher-order 
asymmetries below. A static potential like that illus- 
trated in Fig. |2|c), which is rotated by 45 degrees from 
the RF potentials, is needed to control the relative trap 
strengths in the s and z directions. We add a middle tube 
control electrode between the central cylinder and outer 
tube electrode that we can bias with static potentials to 
create this field. This middle electrode is RF grounded 
with a static potential — Z7 applied to both the top and 
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FIG. 1. (Color Online) (a) A standard four-rod RF linear 
ion trap is transformed into a toroid trap by wrapping the 
electrodes around onto themselves, forming the four rings of 
a toroid trap, (b) An RF toroid trap consists of four ring 
electrodes; RF potentials are applied to the top inner ring 
and the bottom outer ring. This configuration is transformed 
into a compact halo ion trap by contracting the two inner 
rings down to cylindrical electrodes. The outer electrodes are 
extended up and down for convenience in fabrication and as- 
sembly, (c) The resulting form is our new halo trap design. 
Radiofrequency potentials are applied to the top inner elec- 
trode and the bottom outer electrode, creating a circular RF 
node. 



bottom electrodes. Both the inner and outer electrodes 
are also biased with static potentials, -\-U\ and +E/2, re- 
spectively. The combination of the RF pseudopotential 
and this static potential will enable us to vary the trap 
aspect ratio, illustrated in Fig. |2jd). 

There are a large number of free parameters available 
in the design and optimization of this halo trap geometry; 
we optimize a set of these parameters in order to mini- 







FIG. 2. (Color Online) (a) The desired instantaneous RF 
potential is a quadrupole in the s — z plane due to the volt- 
ages zbVb applied to the RF electrodes. Note that the mid- 
dle control electrodes are held at RF ground, (b) The ideal 
pseudopotential is thus symmetric in the s — z plane, (c) 
In order to break the symmetry, static potentials U are ap- 
plied to the both the RF and control electrodes, forming a 
quadrupole rotated by 45 degrees from the RF potential, (d) 
The combination of the static and pseudo potentials gives an 
asymmetric trap, indicating complete control over both the s 
and z trapping potentials. 



mize the higher-order, non-quadratic elements in the RF 
and static potentials. We begin by fixing the inner and 
outer diameters of the cylinder and tubes that make the 
ion trap electrodes. We choose stock parts (from Small 
Parts, Inc.) for these electrodes. The inner cylinder is 
a stainless steel wire with an outer diameter of 510 fim 
(part GWX-0200). The middle tube electrode is a 19 
gauge stainless steel hypodermic round tube with an in- 
ner diameter (ID) of 830 jivn and an outer diameter (OD) 
of 1100 fim (part HTX-19T). The outer tube electrode 
is a 16 gauge stainless steel hypodermic round tube with 
an ID of 1350 /xm and an OD of 1650 (im (part HTX- 
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FIG. 3. (Color Online) Our halo trap design consists of three 
concentric cylindrical conductors. The inner-most conductors 
are angled needle tips, the other two are cylindrical metal 
tubes. The conductors are separated by insulators set back 
from the trap center (not shown). The inner and outer con- 
ductors form the RF trap; the middle conductor allows for 
static voltage control. 



16T). The electrodes are separated by concentric insu- 
lating tubes, matched to the various inner and outer di- 
ameters. We located UHV-compatible stock parts made 
of polyimide that serve as the insulators (parts SWPT- 
028 and TWPT-045). Although we choose specific parts 
for our design and optimization, all of our calculations 
serve for any size trap that has the same ratio of radii 
for the cylinder: (inner tube ID): (inner tube OD): (outer 
tube ID): (outer tube OD) ratio of 1:1.63:2.16:2.65:3.24. 

Having fixed the radii of the three electrodes, the num- 
ber of adjustable free parameters available to optimize 
the potentials is now reduced to the spacing between the 
upper and lower sets of electrodes and the needle tip an- 
gle n (shown in Fig. |2|. We denote the spacings as z n , 
z C: and z t for the inner cylindrical "needle", the middle 
"control" electrode, and outer "tube" electrode, respec- 
tively, as described in Fig. [4] The tip angle will be fabri- 
cated by machining the inner wire prior to trap assembly. 
The other dimensions described in the figure are the in- 
ner cylinder radius R ni the radii of the two tubes R c and 
R t depicted as the distance from the center axis of ro- 
tation to the center of the electrode, and the tube wall 
thicknesses t c and t t for the control and outer tube elec- 
trodes. These five dimensions are fixed by the electrode 
geometry described above. 

We define three dimensionless parameters from these 
physical parameters and vary them to optimize the po- 
tentials. First, the trap aspect ratio Ah is defined as 



the ratio of the average z separation of the needle and 
tube to the tube radius: Ah = (2z n + 2z t )/(2R t ). This 
roughly corresponds to the overall aspect ratio of the 
trap, or the distance between the top and bottom elec- 
trode structures. The second parameter is the "key- 
stone" Kh denned as the ratio of the outer tube sepa- 
ration over the needle separation: Kh = 2z t /2z n . This 
parameter describes how much farther apart the outer 
tubes are compared to the inner cylinder. The final pa- 
rameter is defined as Vh = (z c — z t )/z n and describes 
the control electrode separation as compared to the tube 
separation in units of the needle separation. These pa- 
rameters define several possible configurations simply: if 
all three electrodes are separated by the same amount 
equal to the tube radius R tl the parameter set would 
be: Ah = 2,Kh = 1,V^ = 0. A parameter set of 
Ah = l,Kh = l,Vh = 1 corresponds to separations of 
z n = Rt/2, z c = R t , and z t = Rt/2. The geometry 
illustrated in Fig. [3] corresponds to a parameter set of 
approximately Ah = 0.7, Kh — 1.7, Vh — 2, and a needle 
angle of 17°. 



III. RF AND STATIC POTENTIAL MODEL 

Although motion of ions in an RF Paul trap can be 
described in a variety of ways, the typical approach is 
to model the trapping potentials as an ideal quadrupole 
field [25]. The oscillating trapping potential &(r,z,i) 
for a toroidal trap in cylindrical coordinates r and z 
can be written as a spatially varying potential V(r, z) 
which oscillates at the trap frequency Qt- &( r , z ,t) — 
V(r,z) cos Vt T t. The spatial component of the potential 
is typically described as a quadrupole field, V(r, z) = 



Vo/r% (z 



where Vq is the potential applied sym- 



metrically to hyperbolic electrodes a distance tn from 
the trap center. A toroidal trap like this halo trap 
shifts the trap center radially a distance R from the z 
axis as shown in Fig. |4j Although the purely quadratic 
field described by shifting the quadrupole field V(r, z) = 
V$/r 2 N [z 2 — (r — R) 2 ) is not possible, we are interested a 
trapping potential where the quadratic term is dominant 
[TO] . Thus, we fit the actual halo trap potential to the 
model hyperbolic potential V(r,z) = Vb/^ F [z 2 — s 2 ), 
where s = r — R, and £rf is a single fit parameter. How- 
ever, since the RF fields are rotated by 45 degrees from 
the s — z axes, as seen in Fig. [2ja), we rotate the coordi- 
nate system of the hyperbolic potential in order to com- 
pare it directly with the halo trap field. The quadratic 
potential model describing the field in Fig. [2|a) is thus 



Vrf(s,z) = -j^sz. 



(i) 



The static control potential shown in Fig. [2|c) can also 
be modeled as a hyperbolic potential. Although we ap- 
ply different control potentials to the needle and tubes, 
we fit the actual static potential with a single effective 



(unit) potential C/ e ff and the fit parameter Static- 
corresponding static potential model is then 
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We calculate the actual RF and static potentials, 
VcH(r, z) an d [/ch(^,4 for this compact halo (CH) ge- 
ometry for a given parameter set using finite element 
analysis. We use single-parameter fits (with Vb = 1 V 
and [/"eff = 1 V) to find the normalization distances ^rf 
and Static f° r the RF and static fields. We then define a 
X 2 metric for the two-dimensional scalar fields to deter- 
mine the goodness of the fit. We integrate the square of 
the difference between the real and hyperbolic potentials 
over a circle area A covering the center of the trapping 
region, 

Xrf = // (Vrf (s,z) - Vcn(s,z)f dsdz (3) 

A 

Xstatic = / / (E4tatic(s, z) - U C h(s, z)f dsdz. (4) 



We minimized both of these metrics simultaneously 
through an iterative Monte Carlo approach described be- 
low. 

We also investigate the depth of the RF trap, a partic- 
ular concern for small and novel trap geometries [4j |5] . 
The depth is found by calculating the pseudopotential 
ip(r,z) from the calculated RF trap potentials [26 : 
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where q and m are the ion charge and mass. There is a 
saddle point in the compact halo pseudopotential located 
along the r axis approximately between the outer control 
electrode. The trap depth in eV is shown in Fig. [5] for an 
applied voltage of Vo = 300 V, 24 Mg + ions, and a trap 
frequency of Qt = 27r x 80 MHz. 



IV. HALO TRAP OPTIMIZATION 

We optimize the RF and static potentials using the 
following Monte Carlo procedure: 

1. Optimize the RF potential. 

(a) Choose initial values for the adjustable geo- 
metric parameters Ah, Kh : Vh, an d n . 

(b) Calculate the RF potential using finite ele- 
ment analysis. 

(c) Fit the RF potential to the model and find 

(d) Calculate the Xrf goodness of fit. 
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FIG. 4. (Color Online) There are four key dimensionless 
parameters that define the geometry of the halo trap elec- 
trodes: the needle angle n , the trap aspect ratio Ah = 
(2z n + 2zt)/(2Rt)j the trap keystone Kh = (2zt)/(2z n ), and 
the control electrode offset Vh = (z c — z t )/z n . The other di- 
mensions are fixed and are properties of the conductors and 
insulators used in the model. 
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FIG. 5. (Color Online) The pseudopotential ip for the op- 
timized compact halo trap has a trap depth constrained by 
the saddle point along the z = symmetry line. The pseu- 
dopotential is shown in electron-volts for a 24 Mg + ion with 



an applied potential of Vb 

Qt = 2tt x 80 MHz. 



300 V, and a trap frequency of 



(e) Adjust all of the parameter values and then 
repeat steps (b) through (d) until Xrf * s mm ~ 
imized. 

2. Optimize the static potential. 

(a) Locate the trap center R along the r-axis. 

(b) Calculate the static potential using finite ele- 
ment analysis using the optimized geometric 
parameters. 

(c) Adjust the static potentials /To, E/i, and U2 
until the static potential saddle point lies at 
R along the z = axis. 

(d) Fit the static potential to the model and find 

^static 

(e) Calculate the Xstatic goodness of fit. 

3. Iterate the entire process until both x^ F and Xstatic 
are minimized. 

Since the saddle points for the RF and static poten- 
tials do not necessarily have the same spatial location, 
unlike traditional four-rod traps, we simultaneously op- 
timize both the RF and static potentials under the con- 
straint that the saddle points overlap along the r-axis. 
The optimized RF and static fields are shown in Fig- 
ures [6] and [7| The insets in each figure show a close-up 
of the potentials near the trap center and the residuals 
V RF (s, z) - VchO, z) and U 8ta .tic(s, z) - U C yl(s, z) for the 
RF and static potentials, respectively. The trap center 
is located at R « 430 //m; the normalization coefficients 
are £rf ~ 413 fivn and 4tatic ~ 328 fim. The optimized 
parameters are listed in Table [IJ 



TABLE I. The optimized values for both the geometric pa- 
rameters and the static potentials. 



Parameter 


Value 


Parameter 


Value 


A h 


0.676 


V 


IV 


K h 


1.68 


U 


-42.97 V 


v h 


2.06 


Vx 


1.09 V 


o h 


16.7° 


u 2 


1.03 V 



Both the RF and static potentials have large zones near 
the trap center that approximate very well (to greater 
than 95%) the harmonic potential models. Although it 
may be possible to further optimize the potentials by us- 
ing custom- fabricated electrode structures, the potentials 
generated by stock electrode parts should be of sufficient 
quality for quantum information and Coulomb crystal 
experiments. 



V. CRYSTAL PHASE TRANSITION 

In this final section we derive the first non-trival ion 
crystal phase transition for two ions in the compact halo 
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FIG. 6. (Color Online) The optimized RF instantaneous 
quadrupole potential Vch (with ±1 V applied to the needle 
and outer ring electrodes) across the trapping region. The 
ideal quadrupole fit to the potential has a normalization dis- 
tance of -£rf ~ 413 /mi. The bottom left inset shows the poten- 
tial near the center of the trap with contour lines separated by 
10 mV. The bottom right inset shows the difference between 
V op t and the ideal quadrupole potential Videai with contours 
spaced at 0.5 mV. The central flat region is approximately 
100 /im long by 100 /xm wide, much larger than typical ion 
excursions from the RF node in an ion trap. 



trap. Whereas the first non-trivial phase transition for 
ion crystals in a linear trap occurs with three ions in an 
anisotropic trap [13] , there is a second-order phase transi- 
tion for two ions in the halo trap between the crystal con- 
figuration where both ions lie in the z = plane (Fig. [8]) 
and the regime where both ions are located at stable 
equilibrium positions equidistant from the z = plane 
as shown in the figure. In deriving the condition for the 
phase transition, we follow the derivation of Schweigert, 
et al. [18 with the key difference that we consider a har- 
monic trap in both the r and z diminsions, allowing the 
ions to move in three-dimensional space. 

We consider a three-dimensional system of identically 
charged particles of mass m and charge q confined to 
our halo ring by a ring-shaped external potential V r = 
(l/2)muj 2 .(r — R) 2 (with trap radius R as above) and 
a second, independent, harmonic potential in the z- 
direction V z = (l/2)muj 2 z 2 • The interaction Hamilto- 
nian, including Coulomb repulsion, of this classical sys- 
tem of N particles is thus 
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FIG. 7. (Color Online) The optimized static quadrupole po- 
tential Uch across the trapping region. The ideal quadrupole 
fit to the potential has a normalization distance of Static ~ 
328 /mi. The bottom left inset shows the potential near the 
center of the trap with contour lines separated by 2 mV. The 
bottom right inset shows the difference between Ucu and 
the ideal quadrupole potential [/static with contours spaced 
at 0.5 mV. The central flat region is approximately 100 /mi 
long by 100 /im wide, much larger than typical ion excursions 
from the RF node in an ion trap. 
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FIG. 8. Two ions in the halo trap begin on opposite sides of 
the ring. There is a second-order phase transition where the 
ions shift off of the ring in opposite directions. 



for each particle located at r^. Following the convention 
in [18] , we re-scale the length and energy units in order to 
work with dimension- less parameters. The length scale 
is denned as the ratio of the Coulomb and radial trap 
coupling constants, 



( q 2 2 
\47reo muj% 
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The length scales for a number of different common ion 
trap species are listed in Table \H\ All of these are similar 
in scale the model compact halo trap radius. We also 
define an energy scale based on the radial potential and 
this length scale, E* = (l/2)muj^r^ and scale the overall 
energy by this parameter. Finally, we define the trap 
aspect ratio as the ratio of the z and r trap frequencies, 
similar to a linear trap, 



ut 
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TABLE II. Various scaled radii for common ion trap config- 
urations. The final entry is for 300 nm diameter polystyrene 
spheres, charged to 1.6 x 10 -16 C. 



Ion 


UJr/2-K 


r* 


24 Mg+ 


2 kHz 


419 fim 


4u Ca+ 


1.5 kHz 


428 /mi 


m Yb + 


0.8 kHz 


401 /xm 


300 nm PS 


100 Hz 


429 /mi 



In the simplest configuration, the ions all lie in the 
z = plane and distribute themselves equally around a 
circle [12 with their center of mass located at the center 
of the ring x crn = y cm = z crn = 0. Since no external 
forces act on the ions in this trap, the center-of-mass of 
the system will not move. We use this to simplify the 
ion locations for the N = 2 case, making the assumption 
that the ions start on opposite ends of the x-axis, the 
positions of the two ions will be related by the following 
three constraints: 



X\ = — X2 
2/1=2/2 = 
z\ = -z 2 . 

These assumptions reduce Eq. (|6| to 

1 



H 



2V^ 



2az 2 + 2(\x\ -r ) 2 . 



(9) 



(10) 



where we have scaled the distances, x = x\jr^ z = z\jz±, 
and energy H — >• H/E* as noted above. The trap center 
is also scaled by the same parameter, r*o = R/r*. This 
unitless radius now describes the ratio of the physical 
trap R to the strength of the Coloumb interaction. We 
will assume that the trap is small enough that the ions 
interact strongly across the diameter of the trap circle. 

The equilibrium positions of the two ions are found 
from the first spatial derivatives of the interaction Hamil- 
tonian in the x and z directions, 



dx 
dz 



=0 



=0. 



(11) 
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FIG. 9. (Color Online) (a) The scaled z position as a function 
of a and tq. The phase transition occurs along the line de- 
scribed by Eq. (13). (b) The scaled x position also shows the 
phase transition, (c) The phase transition in a and ro space. 



The first set of stable equilibrium positions are at 



16 l/3 r 2 



3 (27 + 16rg - 3^81 + 96^) 
(27 + 16rg - 3V81 + 96rg) 



1/3 



432V3 



=0. 



(12) 



This solution is independent of the trap aspect ratio a 
and is valid for 



\ a ~ l \ A 

r o > 1/0 and a > 1. 



(13) 



As the ring becomes small, or the trap aspect ratio de- 
creases, corresponding to a weakening of the z potential, 
the ions reach a point where they shift off of the z = 
plane into the second crystal phase described by the equi- 
librium positions 
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1 /(a-l) 2 -4r 2 a 2 / 3 



(a - l) 2 a 2 / 3 



(14) 



The x and z positions as a function of ro and a are shown 
in Fig. [9ja-b). This is a second-order phase transition 
since the derivative of the mean radial position of the ion 
cloud is discontinuous. The boundary between the two 
phases is also shown in Fig. |9jc). 

We have shown that there is a non-trivial phase tran- 
sition for two ions in the compact halo trap. Future work 
includes solving for the phase transition for larger num- 
bers of ions, in particular for odd numbers of ions, which 
have a different behavior for small numbers of ions [18] . 
We have also presented an optimized electrode configu- 
ration for creating the compact halo trap using readily- 
available electrode parts. This type of trap is an alterna- 
tive topology to more common point and linear ion traps 
and may be useful as a small, novel trap geometry for 
some quantum information applications. 
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